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Abstract. The powerful wind- wind collision in massive star binaries creates 
a region of high temperature plasma and accelerates particles to relativistic en- 
ergies. I briefly summarize the hydrodynamics of the wind-wind interaction and 
the observational evidence, including recent 7-ray detections, of non-thermal 
emission from such systems. I then discuss existing models of the non-thermal 
emission and their application to date, before concluding with some future 
prospects. 



1 Introduction 

The winds of O and WR stars couple high massdoss rates of ~ 10 -7 to a few 
times 10~ 5 MQyr" 1 with fast velocities typically ~ 1000 — 3000 kms" 1 . The 
ram pressure balance between the hypersonic winds determines the position of 
the wind-wind collision region (WCR). Shocks are formed either side of the 
WCR which thermalize the plasma, heating it to ~ 10 7 — 10 8 K. The proper- 
tie s (e.g. density) of the WCR can span a very large range (see e.g. Table 2 
in 



Pittard et al. 20051 ). reflecting the diversity of the underlying binary popu- 



lation. In systems where the orbital period is only a few days, the shocks are 
collisional, and the WCR displays a large aberration and downstream curvature 
due to the coriolis force. Cooling of the shocked gas is also likely to be signifi- 
cant. In addition, the presence of the companion star's radiation field may have 
significant influence on the driving and dynamics of the winds. In contrast, in 
systems with longer orbital periods the shocks may be collisionless, orbital and 
radiation-field-induced effects on the dynamics and geometry of the WCR are 
much smaller, and the hot shocked gas may only cool adiabatically as it flows 
downstrea m. Particle acceleration at t he shocks (but also possibly within the 
WCR - see lPittard k, Dougherty! ( 20061 ) for a detailed review) is likely to occur 



in these wider systems. Systems with (highly) eccentric orbits are particularly 
interesting since the changing separation of the stars is useful as a probe of the 
physics which takes place. Colliding wind binaries (CWBs) are also useful as 
simpler, less complicated, analogues of the multiple wind-wind collisions which 
occur throughout the volume of clusters of massive stars. 

The collision of t he winds is best studied using hydro dynamical simulations. 
A seminal paper by IStevens. Blondin fc Poltockr (119921 ) was the first to self- 



consistently include cooling, and it focussed on the dynamical instabilities and X 
ray emission which arise from the WCR. Other works since then have examined 



the influence of th e companion star's radiation field (IStevens fc Pollockl [ 1994; 



Gavlev et al. 19971 ). the effects of thermal conduction ( Mvasnikov fc Zhekov 
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19981 ). electron-ion tempera ture non-equi librium ( Zhekov fc Skinner 20001 ). and 



non-equilibrium ionization ( Zhekov 20071 ) . The effect of clumps within the stel- 
lar winds on a WCR ne ar th e radiative and adiabatic limits has been inves- 
tigated by IWalderl ( 1998) and iPittardl ( 20071 ). resp ectively. Three-dime nsional 
models have allowed orbital effects to be explored (jLemaster et al.ll2007l ). This 
work has also been comple mented by smoothed-particle-hydrodynamics simula- 
tions ( Okazaki et al. 20081) and the development of a co mpute-efficient dynam- 
ical model ([Parkin fc Pittardl 120081 : iParkin et aT1l2009l ). More recently, three 
dimensional hydrodynamical models which incorporate gravity, the driving of 
the winds, the cooling of the shocked plasm a, and the orbital motion of the 
stars, have been presented by IPittard ( 20091 ). This work represents the most 
realistic simulation to date of the wind-wind collision in short period systems. 

Fig. Q] displays cuts through the orbital plane of the density and temper- 
ature in two different systems with identical O-stars but with orbital periods 
of 3 (left) and 10 days (right). The WCR is highly radiative in the shorter 
period system since the winds collide before they have had much opportunity 
to accelerate, the postshock gas temperature is lower, and the postshock den- 
sity is higher. There is also considerable inhibition of the acceleration of the 
winds caused by the opposing radiation field of the companion star. In contrast, 
with an orbital period of 10 days the wider separation of the stars is enough to 
more than double the pre-shock velocities of the winds (along the line-of-centres 
through the stars it increases from 730 kms" 1 to 1630 kms^ 1 ), which leads 
to significantly higher (lower) postshock temperatures (densities), resulting in 
a largely adiabatic WCR. The higher wind speeds and reduced orbital veloci- 
ties also diminishes the aberration and downstream curvature of the WCR due 
to coriolis forces (note the difference in spatial scale between the two systems 
presented in Fig. [T]). 



2 Particle Acceleration in Colliding Wind Binaries 

Direct evidence for no n-thermal emis s ion fr om relativistic particles within a 
WCR was presented bv lWilliams et a~ ( 19971 ). Here, radio images were overlaid 



on UKIRT shift-and-add IR images of WR 147. When the southern (thermal) ra- 
dio source was aligned with the southern (WR) star, the northern (non-thermal) 
radio source was found to lie just south of the northern (O) star, in a position 
consistent with the point of ram-pressure balance between the win Direct 
imaging of the WCR i n WR 146 and WR 140 have provided further support 
for th is interpretation ( Dougherty. Williams Pollaccol [2000l ; [Dougherty et al 



20051 ). See Dougherty (these proceedings) for further details. 

WR 146 and WR 147 are both very wide, and the non-thermal radio emis- 
sion from the WCR escapes easily from the system. This contrasts with WR 140 
(the best studied of any WR+O binary), where the stars are very much closer 
together, and which displays dramatic, phase-repeatabl e, variations in its ra- 
dio emission over the course of every 7.94 yr orbit (see IWhite &: Becker 19951 ; 



WR 147 is the o nly CWB to date which also has spatially resolved X-ray emission - see 
IPittard et all (|2002l ). 
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Figure 1. Density and temperature plots through the orbital plane of hy- 
drodynamic models of the wind-wind collision in short period 06V+06V 
binaries. In the left plots the orbital period is 3 days, and the winds collide 
at relatively low speeds. Hence, the WCR is highly radiative and is strongly 
distorted by the orbital motion of the stars. In the right plots the orbital 
period is 10 days, the winds collide at higher speeds, and the WCR is largely 
adiabatic. The top panels display the density and the bottom panels the tem- 
perature. The left panels have si des of length 24 R Q , while the right panels 
have sides of length 570 R Q . See lPittardl (|2009t) for further details. 



4 



J. M. Pittard 



Dougherty et al 



20051 ). The high eccentricity of the orbit causes the stellar 
separation to vary between ~ 2 — 28 AU. While the radio lightcurve has defied 
satisfactory explanation to date, it is likely that at least part of its variation 
is caused by variable circumstellar extinction to the source of the non-thermal 
emission as the O star orbits in and out of the radio photosphere in the dense 
WR wind. More recently, this system has been imaged with the VLBA, yielding 
a full orbit definition, i ncluding, most importantly, the inclination of the system 
( Dougherty et alj|2005l ). While this study has helped to provide some of the best 
modelling constraints of any system, relatively little is known about the wind of 
the O star, and the wind momentum ratio of the system remains ill-constrained. 

In the hard X-ray and 7-ray regimes there has been many false dawns con- 
cerning potential emission from CWBs. However, definitive evidence for non- 
thermal X-ray and 7-ray emission from CWBs has finally been presented in 
recent years, wit h an especially exciting p e riod i n just the last few months. Us- 
ing INTEGRAL, iLevder. Walter fc Rauwl (|2008l ) presented a clear detection of 
MeV 7-ray emission from 77 Carinae (a lo wer spatial resolut ion observation with 
BeppoSAX was previously presented by IViotti et al. 2004 , t h ough some of the 
emission attributed to r/ Carinae was shown by Levder et al. ( 20081 ) to be asso- 
ciated with other near by sources). This dete ction has since been confirmed with 
a Suzaku observation ( Sekiguchi et al.l [20091 ). And i n just the last few months, 
77 Carinae has been associated with an AGILE source ( Tavani et al. 20091 ). and is 
on the FERMI bright source list. Upper limits h ave also been placed on the MeV 
emission from WR140, WR146, and WR14 7 dDe Becker et aTl l2007f). and on 
the TeV emission from WR 146 and WR 147 (|Aliu et al.ll2008l ). Finally, a great 
deal of effort has gone into searching for non-t hermal X-ray emi s sion from CWBs 
with known non-thermal radio emission (e.g. Rauw et al. 20021 ; |Pe Becker et al ' 



20041 . 12006ft . However, the lack of success in this undertaking appears to show 



that the physical conditions necessary to produce observable emission at radio 
and hard X-ray energies are mutually incompatible. 

There are many possible mechanisms for accelerating particles in CWBs. 
Most works assume that diffusive shock acceleration (DSA) is the dominant 
process, though reconnection and turbulent processes are other possibilities. 
Each mechanism differs in its efficiency, and in the properties of the non-thermal 
particles which it produces, such as thei r spectral index. A de t ailed discussion 
of the many possibilities can be found in lPittard Dougherty! ( 20061 ) . 



3 Models of the non-thermal emission 



3.1 Early models 



So me notes on the key p hysics of particle acceleration in CWBs were presented 
bv lEichler fc Usov ( 19931 ). Early models of the non-thermal radio emission were 



very simple, with the observed flux (S° hs ) assumed to be a combination of the 
free- free flux from the spherically symmetric winds (S^), plus the flux from 
a point-like non-thermal source located at the stagnation point of the winds 
(S^), the latter being attenuated by free-free absorption (opacity r^) through 
the surrounding winds: 

S?» = £fi + S?e- 1 ?. (1) 
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This approach allows relatively si mple analytical solutions to the radiative trans - 
fer equation to be obtained (e.g. IWilliams et al.l [l99Cll ; IChapman et al.l Il999l ). 



Unsurprisingly, such simple models fail t o reproduce the spect ral variation of 
the emission with orbital phase. This led Williams et al. ( 19901 ) to propose the 



need for a more complicated model which accounted for the lo w-opacity "hole" 
i n the dense WR wind created by the O-star's wind. However, IWhite &: Becker 



( 19951 ) pointed out that in the case of WR140, even the O-star's wind has sig- 



nificant opacity. These works clearly demonstrate the need for more realistic 
modelling where the presence of the hot, low-opacity, WCR is a fundamental 
component. In addition, the interferometric observations of spatially-extended 
synchrotron emission indicate that the assumption of a point-like non-thermal 
source needs also to be modified. Other problems with these simple models are 
that no consideration of cooling (e.g. inverse Compton) or other mechanisms 
(e.g. the Razin effect) are made. 

In many ways, the state of modelling in the 7-ray domain was, until a few 
years ago, even more rudimentary. With no firm detections at 7-ray energies, 
theoreticians were reduced to estimating the inverse Compton (IC) luminosity, 
Li c , by the following simple formula: 



r. - E^Lt (o\ 

^ic — t j J -'sync> 



where L sync is the synchrotron luminosity, and U p h and Ub are the photon and 
magnetic field energy densities, respectively. A fundamental problem with the 
use of the above formula is that the predicted value of L K is highly sensitive to 
the assumed magnetic field (B = \ Z8irUn). Varying B r esults in a wide range 
of predictions for L 1C , as shown in iBenaglia fc Romero! ( 20031 ). The magnetic 



field in the WCR is highly uncertain, even when the surface magnetic fields are 
known (which is very rare), because one cannot necessarily extrapolate these to 
obtain their strength in the WCR, since there are the possibilities of large-scale 
magnetic reconnection within the WCR on the one hand, and magnetic field 
amplification by non-linear DSA on the other. Finally, observed rather than 
intrinsic values have been used for L syQC . This may be of little consequence in 
the wider systems where the attenuation of non-thermal emission through the 
circumstellar envelope surrounding the stars may be negligible, but for a given 
set of parameters, L lc could be underestimated in closer systems. 

While there are many problems and uncertainties associated with the early 
modelling work of the non-thermal radio and 7-ray emission from CWBs, these 
simple models nevertheless paved the way for the more sophisticated modelling 
which has since followed, as we now describe. 

3.2 Recent developments 

3.2.1 Models of the radio emission 

A major step along the path toward s improved model s of th e radio emission from 
colliding wind binaries was taken bv lDoughertv et al. ( 20031 ). This work removed 



the assumptions of a point-like source of non-thermal emission, and a spheri- 
cally symmetric, single temperature, surrounding envelope. Instead, models of 
the thermal and non-thermal radio emission were based on 2D, axisymmetric hy- 
drodynamical simulations. This approach provided a much better description of 
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the density and temperature distribution in the the system, allowing sight-lines 
to the observer to pass through regions of both high and low opacity. The non- 
thermal emission was treated in a phenomenological way: accelerated electrons 
were assumed to be present within the WCR, with an energy density (U ie i >e ) 
proportional to the local thermal energy density (U t h) i.e. J7 r ei,e — Crei.e^th- The 
magnetic field energy density was specified in a similar manner: U-q = CB^th- 
The non-thermal electrons were further assumed to have a power-law distribu- 
tion, N(j)dj = C^~ p d'Y, where 7 is the Lorentz factor and C is proportional to 
Crei,e- The non-thermal electrons were assumed to arise from test particle DSA, 
where p = 2 for strong shocks of adiabatic index 5/3. 

Though this work introduced a number of assumptions and involved some 
essentially free parameters, it nevertheless provided a great deal of new insight 
into the phenomenon of radio emission from colliding wind binaries. An imme- 
diate benefit was the realization of the potential importance of the Razin effect 
in attenuating the low frequency synchrotron emission within the WCR. It also 
established several key scaling relations. For instance, given the assumptions 
in the model, the total synchrotron emission from the entire WCR in adiabatic 
systems was found to scale as D -1 / 2 ^" 1 / 2 , where D is the separation of the 
stars (for comparison, the X-ray emission in the optically thin, adiabatic limit, 
scales as D~ l ). This work also highlighted the importance of IC cooling, which 
was noted to be important even in wide systems. Indeed, neglect of IC cooling 
of the non-thermal electrons in this model led to an overestimation of the high 
frequency synchrotron flux in a model of WR 147. 

This failing was addressed in a follow-up paper (jPittard et alJl200d ). Here, 
the non-thermal electrons were assumed to be accelerated by DSA at the global 
shocks confining the WCR, and to cool once they flowed downstream. Since 
the amount of cooling is dependent on the "exposure" time of the non-thermal 
electrons to the radiation fields of the stars, this assumption led to those elec- 
trons suffering the most severe cooling to concentrate along the contact surface 
between the winds, and thus to a dearth of emission from this region (see Fig. [2]). 
This addition significantly improved the fit between models and observations of 
WR 147. It also modified the scalin g relation for the total synchrotron luminos- 
ity noted in Dougherty et all ( 20031 ) - now, the intrinsic luminosity was observed 
to decline with stellar separation as IC cooling became increasingly strong. In 
addition, it was noted that the thermal radio emission from the WCR scales as 
-D -1 , in an identical way to the thermal X-ray emission. Since this emission is 
optically thin (on account of the high temperatures within the WCR), it can 
mimic a synchrotron component. Therefore, one needs to cautiously interpret 
data with a spectral index —0.1 < a < 0.5, where S v oc v a (more negative values 
of a, e.g. a ~ —0.5, clearly indicate a bona-fide synchrotron component). 



3.2.2 Models of the non-thermal X-ray and 7-ray emission 

The dramatic sensitivity gains made by arrays of Cerenkov telescopes have led 
in recent years to a new interest in t he level of non-t hermal X-ray and 7-ray 
emission from colliding wind binaries. iBednarekl ( 20051 ) calculated the expected 
7-ray emission from WR 20a, a WR+WR binary. The short orbital period means 
that the optical depth to electron-positron pair creation is high enough to initi- 
ate electromagnetic cascades in this system. Particle acceleration by magnetic 



Non-thermal emission models of early-type binaries 



7 




Figure 2. The impact of IC cooling on the intensity distribution of a CWB 
model fsee lPittard et al1l2006l for details), for a viewing angle of 0° and at 
1.6 GHz (top) and 22 GHz (bottom). The images on the left do not include 
IC cooling, while those on the right do. Each image has the same intensity 
scale and contours. 
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Figure 3. a) Geometry of the 2-zone model in iReimer et alJ (j2006(). b) 
Evolution of the non-thermal electron spectrum from the inner acceleration 
zone (solid) line as a function of downstream distance in the advcction zone. 
At low energies adiabatic/expansion losses dominate, while at high energies 
IC losses dominate, c) As b) but for nucleons. Only adiabatic losses occur. 
See IReimer et al.l (|2006f) for further details. 



reconnection and DSA was considered, and it was concluded that detectable 
neutrino fluxes should be produced. Due to the high optical depth to TeV pho- 
tons with in this system, WR20a cannot be directly responsible for nearby TeV 
emission ( Aharonian et af]|2007l ). which more likely is the result of acceleration 
processes within the collective win d of the nearby cluster Westerl und 2. 

Following on from this work, Reimer. Pohl fe Reimer ( 20061 ) developed a 
two-zone model of the non-thermal emission. Particles are accelerated in an inner 
zone where their spatial diffusion exceeds their motion due to advection with 
the background fluid. Their energy distribution is self-consistently computed 
by considering all relevant gain and loss mechanisms. Particles are assumed to 
be resident within this region until their timescales for advection and diffusion 
are comparable, after which they are assumed to move into the advection region 
where they suffer further losses as they flow downstream. Fig. [3] shows the 
assumed geometry and the resultant non-thermal energy spectra of the electrons 
and nucleons. 

Consideration is also given to the anisotropic nature of the IC process, 
where the emitted power is dependent on the scattering angle. Fi g. [H shows the 



predic ted IC emission from WR 140 as a function of orbital phase. IReimer et al 
(120061 ) conclude that while WR 140 should be easily detected with GL AST/Fermi 
the change in the IC flux with viewing angle due to anisotropic scattering is likely 
to be obscured by large variations in the energy density of the stellar radiation 
fields resulting from the high orbital eccentricity. However, the latest work by 



Reimer fc Reimerl ( 20091 ) demonstrates that it is possible to use the property of 



nonisotropic IC emission to constrain the orbital inclination of colliding wind 
systems. 

Another model of the X-ray and 7-ra y emission, which i s in many ways com- 
pleme ntary to that of IReimer et al.l ( 20061 ). was presented bv lPittard fc Dougherty! 
( 20061). This bu i lt on the phenomenologi c al mo del developed previously by 
iDoughertv et alJ (120031 ) and iPittard et at! (|2006l ) to explore the non-thermal 
radio emission. Although the energy spectrum is as sumed rather than calcu- 
lated, and so in this sense is weaker than the model in IReimer et alJ (120061 ). this 
approach benefits from a realistic description of the density and temperature dis- 
tribution within the system, and constraints placed on the key parameters (e.g. 
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Figure 4. Predicted IC spec tra for WR140 at phases 0.2, 0.671, 0.8 and 
0.955 from lReimer et al.l (|2006f ). 7-ray absorption is not included. 



mass-loss rates) by fits made to the X-ray data. IPittard fc Dougherty! (|200fll ) 
showed that other uncertainities, such as the particle acceleration efficiency and 
the spectral index of their energy distribution (both of which unfortunately re- 
main ill-constrained from fits to radio data), have at least as much influence on 
the predicted flux as the angle-dependence of the IC emission. 

Fig. [5] shows a predicte d spectral energy distr i butio n for WR 140 from one 
of the models presented in IPittard fc Dougherty! ( 20061 ). Large differences in 
the predicted 7-ray emission occur depending on whether the low frequency 
turndown in the radio spectrum results from free-free absorption through the 
surrounding stellar winds, or from the Razin effect. Furthermore, satisfactory 
fits to the radio spectrum at orbital phase 0.837 could be achieved in one of two 
ways: either with a standard p = 2 index, or with a harder index (e.g. p ~ 
1.4). Harder indices can result from the shock re-acceleration process, whereby 
the n on-thermal particles pass through a sequence of shocks ( Pope &: Melrosel 
Il994j ) . or from 2nd order Fermi acceleration. Either of these processes may 
be significant in CWBs, since the clumpy nature of the winds means that the 
W CR is likely to be highly turbulent, with weak shocks distributed throughout 



it (Pittard 



While 



20071) . 



Pittard fe Dougherty! ( 20061 ) could not determine the nature of the 



absorption process from the fits to the radio spectrum, future 7-ray detections 
will determine the 7-ray flux and spectral index, and thus will also distinguish 
the nature of the low-frequency turndown. The acceleration efficiency of the 
non-thermal electrons and the strength of the magnetic field will then both be 
revealed. 
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Figure 5. The radio a nd non-thermal UV, X-ray and 7-ray emission cal- 
culated from model B in iPittard fc Doughertvl (|2006T ) , together with the ob- 
served radio and X-ray flux (both at orbital phase 0.837). The model IC (long 
dash), relativistic bremsstrahlung (short dash), and neutral pion decay (dot- 
ted) emis^ion_romp_caients are s hown, along with the total emission (solid). 
See IPittard fc Doughertvl (|2006t ) for more details. 



Work is also continuing on fits to the radio spectra of WR 140 at other 
phases around its orbit. Preliminary results indicate that there is significant 
evolution of key parameters in the model (such as the acceleration efficiency 
and magnetic field). F inally, the non -thermal emission from the short period 
O+O binary models in Pittard ( 20091 ) is being investigated. One area of inter- 
est concerns the ease with which non-thermal radio emission can escape these 
systems, despite the orbital-induced curvature preventing lines-of-sight to the 
non-thermal emission from existing purely within the low opacity, hot WCR. 



4 Conclusions 

In recent years there has been a steady improvement in theoretical model pre- 
dictions of the non-thermal emission from CWBs, in both the radio and 7-ray 
domains. This work is being driven by corresponding advances on the obser- 
vational front. In the past, CWBs have played a poor second role to SNRs in 
terms of investigations of the physics of high Mach number shocks and cosmic 
ray acceleration. However, they provide access to higher mass densities, radia- 
tion backgrounds, and magnetic field energy densities than found for supernova 
remnants, thus enabling studies in previously unexploited regimes, and due to 
recent observational detections of non-thermal emission at hard X-ray and 7- 
ray energies are slowly gaining popularity in the community. Future high energy 
observational prospects with Fermi and CTA, etc., look very good. An exciting 
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future is also expected in the radio regime, with EVLA, e-Merlin, and SKA (see 
Benaglia, these proceedings, for further information). 

Future theoretical work should combine the best features of the different 
modelling efforts to date: for instance, adding a calculation of the non-thermal 
pa rticle energy spectra and ani sotropic IC to the hydrodynamical based models 
of lPittard & Dougherty! <|200d ). A more distant goal would be to incorporate 



the effects of particle acceleration on the underlying thermal plasma, since DSA, 
when present, appears to be very efficient at placing energy into the non-thermal 
particle distribution. 
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